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Abstract 

We consider the canonical symplectic form for sine-Gordon evaluated explicitly on the solitons of the 
model. The integral over space in the form, which arises because the canonical argument uses the 
Lagrangian density, is done explicitly in terms of functions arising in the group doublecrossproduct 
formulation of the inverse scattering procedure, and we are left with a simple expression given by 
two boundary terms. The expression is then evaluated explicitly in terms of the changes in the 
positions and momenta of the solitons, and we find agreement with a result of Babelon and Bernard 
who have evaluated the form using a different argument, where it is diagonal in terms of 'in' or 'out' 
co-ordinates. Using the result, we also investigate the higher conserved charges within the inverse 
scattering framework, check that they Poisson commute and evaluate them on the soliton solutions. 
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1 Introduction 



The purpose of this paper is to calculate the symplectic form restricted to the phase space describing 
the dynamics of solitons in the sine-Gordon model, in terms of the loop groups of the inverse scattering 
procedure. 

The aim is to be able to do Hamiltonian mechanics with solitons, while retaining all their positions 
and phases. The reader may remember that any method proceeding only via the monodromy matrix, 
such as the quantum inverse scattering method looses information on the positions and phases. 
Some work has already been done on the soliton symplectic form restricted to the solitons of sine- 
Gordon by Babelon and Bernard jjj, where a result was obtained in terms of the positions and 
momenta of the solitons by using an intuitive argument based on 'in' and 'out' co-ordinates; the form 
becoming diagonal in these co-ordinates. This argument solved for the form without evaluating it 
explicitly (the integrals over space were not performed). 

The approach adopted in this paper relies on an interesting simplification which occurs when the 
symplectic form is expressed in the language of inverse scattering and group factorisations. The space 
integral in the form can be done explicitly by an 'integration by parts trick', and the form is then a 
difference of two boundary terms. This general expression for the form ( 5.1 3| ) is of interest by itself, 



and could be useful for applications such as quantisation - in fact quantisation is the real motivating 
force behind our construction. However, we go further and calculate it explicitly in terms of the 
changes in position and momenta, where we indeed find agreement with the results of Babelon and 
Bernard Q. 

The second purpose of the paper is to relate this to an abstract integrable system consisting of a 
group factorisation and a classical vacuum map. From this abstract point of view we do Hamiltonian 
mechanics by studying Hamiltonian functions arising as 'higher momenta' in integrable theories, and 
we calculate the momenta for the sine-Gordon case. We explicitly verify a result of Olive and Turok 
||, restricted to sine-Gordon, that the values of the higher charges with odd Lorentz spins are zero. 

2 Preliminaries 

We firstly define light cone coordinates on R > by x± = t ± x, and let d± denote differentiation with 
respect to x± respectively. In these coordinates the sine-Gordon equation is written as 
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d+d-u= — — sin(/3u), (2.1) 
P 

where u(x+,X-) is a real valued function on R ' . 

Now consider the following linear system for a function 

* : R 1 ' 1 x C* — > GL 2 (C) 

d + V(x + ,x^, A) = VA(x + ,X-,\), d^(x + ,x^,\) = VB(x + , X-,\). (2.2) 
The linear system is overdetermined, and the compatibility condition derived from 



is precisely 

We shall take A and B of the form 



d-A-d+B =[A,B}. (2.3) 



(3u flu 
A = -(3(d+u)s 3 + 2mA(cos( — )si - sin( — )s 2 ) 

= -/3{d + u)s 3 + 2m\e-' 3uS3 s 1 e t3uS3 
B = I3(d-u)s 3 + 2mA- 1 (cos(^) Sl +sin(^) S2 ) 



1 



= P(d-u)s 3 + 2m\~ 1 e PuS3 s 1 e-P US3 , (2.4) 

where A can be an arbitrary complex number and is called the spectral parameter, and s±, S2, and S3 
are the standard anti-hermitian basis elements of the Lie algebra sl(2) 

[si,s 2 ] = ^s 3 , [s2,S3] = ^si, [s 3 ,si] = is 2 - (2.5) 

Observe that the spectral parameter A is in the principal gradation of the underlying loop algebra. 
The compatibility conditions (2.3) then immediately yield 
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TYl 

{d + d_u H — — sin(/3u))s 3 = 0, 

P 



in other words u is a solution to the sine-Gordon equation (2.1). 
For convenience we shall want to have S3 diagonal, so we choose the explicit basis 



4 v* oy ' " 4 \ -i j ' 4 v° -i 

We have a simple vacuum solution to these equations given by u constant with cos(-|) = 1. In this 
case we define J and K by 

A = J = 2m\s 1 , and B = K = 2m\' x s 1 . (2.7) 

The equations 

9 + *o = ^qJ, and d-^ = ^ K 

have a simple exponential solution, and if we then subtract off this vacuum solution from the general 
solution ^ by defining 

= ^o^, (2-8) 

we see that obeys the equations 

<9+0 = 4>A- J0, and <9_0 = 4>B - Kef). (2.9) 
We demand that A = and A = 00 are regular points of the function 0(A), by requiring that 
(jr 1 d+(j) = A-0 _1 J0, and 0~ 1 cL0 = B - 0" 1 K<j) 
are free from singularities at and 00, and from this we deduce that 

e -(h»3 aie pua a = ^(oq) -i Sl 0(oo) ; (2.10) 

and 

e /3us3 Sie -/3« S3 = ^-^^(O). (2.11) 



We can solve (2.10) by setting the normalisation of 0(A), that is 0(oo) = e^ US3 . We can then compute 
the explicit soliton solutions from 0(0) and Q2.1ip , which must satisfy 0(0) = Ne~P US3 , for N a 
constant matrix such that N^^-siN = si. 

To be consistent with the linear system we impose the condition, as shown in H, that 

Cty(A)tft = /(A)0(-A), (2.12) 

where /(A) is a scalar function and U = 4s3. 

It is well known, for example see j^, |7|, or |2| [|, that the meromorphic solutions 0(A), which are 
unitary on the real axis, can be written as products of 

(V + ^^l (2.13) 
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where i runs from one to the number of solitons in the system. Here Pi is a Hermitian projection, i.e. 
Pf = Pi, and Pj = Pi. The product can be taken in any order, but then different P's will apply. The 
positions of the solitons are encoded into the Pi's, and the momenta of the i th soliton is related in a 
simple way to |aj|, that is the rapidity, 9i, is given by |a,| = e~ 6t . 

The condition ( 2.1 2| ) restricts the positions of the poles cti to either the imaginary axis, which 
corresponds to solitons and anti-solitons (corresponding to a single factor of the form ( |2.13 )), or if not 
on the imaginary axis, the poles come in pairs, one at oii and the other at —on. This latter situation 
corresponds to a breather - a bound state of a soliton and anti-soliton. 



3 The doublecross product method. 

At this point we must beg the readers indulgence as we consider a very abstract picture. The justifi- 
cations for our digression are: 

(1) We will derive concrete formulae which we will apply to our present case, the sine- Gordon 
model. 

(2) We shall find general formulae for the 'higher conserved momenta' typical of 1 + 1 dimensional 
integrable field theories. 

(3) The expression found for the symplcctic form will be based on functions defined in this section. 
Consider a group factorisation X = QAi = A4G, where Q and Ai are subgroups of X, and Q n M 

consists only of the identity element. We require that any element of X can be uniquely factorised 
as cup, where c € Ai and ip £ Q, and can be uniquely factorised as (f>f, where / e M and <j> G Q. 
This is the definition of a group doublecross product, and appears in the construction of Hopf algebra 
bicrossproducts however we shall not consider Hopf algebras in this paper. 

To construct an example of an abstract integrable field theory, take a spacetime S, and a function 
a : S — > Ai, which we shall call the classical 'vacuum' map. The group Q may be called the classical 
phase space. For a given classical solution <f>o in the phase space, we perform a factorisation for all 
space-time positions seS; 

a(s) = 4>{s) b(s) , 0(a) e Q , b(s) G Ai . (3.1) 

The solution to the field theory the fields at the point s € S, is encoded into b(s) and </>(s). If S is a 
differential manifold we can recover a linear system from the factorisation ( |3.1| ) by differentiating it 
along a vector g in S to give a g <po = <p g b + <f>b g (we use subscripts for differentiation). This can be 
rearranged to give 

4> g = a g a~ l <j) — ipbgb^ 1 . (3-2) 

For the sine-Gordon model, we take S to be the ordinary 1 + 1 dimensional flat Lorentzian space- 
time. The group Ai consists of functions from C* to GL2(C) which are complex analytic (so that 
essential singularities are likely to be present in elements of Ai at and oo), and are unitary on R*. 
The group Q consists of meromorphic functions from Coo to GL/2( < C) which are unitary on Moo, and 
which satisfy the symmetry condition ( [2.12 ). The solution 0(A) to the linear syst em (2.S) is valued in 



Q, and is the same (j)(X) which appears in the factorisation (3.1). Equation (3.2) has the same form 



as (2.9), provided that a + a = —J, a_a = —K, b + b = —A, and 6_6 = —B. Therefore the 



'vacuum' map o(A, x, t), in this case, is defined as 

a(X,x,t) = e - Jx +- Kx ~ . (3.3) 

The proof of the corresponding statement for the A n _i affine Toda field theories can be found in 
0. Sine-Gordon can be recovered from the analysis in if we restrict ourselves to the case n = 2. 
In the more general affine Toda case the unitarity condition has to be dropped, and a generalisation of 



(2.12) applies. The loss of unitarity means that the residue Pi of the poles in each of the simple pole 
factors ( p. 13 ) is no longer Hermitian. For sine-Gordon the specification of the kernel of Pi is fixed in 



terms of the image, therefore the problems pointed out in H d ue to the kernel do not apply. 
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4 Explicit calculation of some soliton solutions 



The explicit form of the soliton solutions are extremely familiar, and have been known for quite some 
time, however we shall need to know precisely how the positions of the solitons are encoded into the 
residues Pi and for this we need to calculate some solutions using the method. We will also use the 
group factorisation (|3.1|) to our advantage, to simplify the calculations. 



We write 4>(X) as a product of meromorphic loops, of the form (2.13), which we call ip(X), times 
the non-trivial normalization matrix 

g = <l>(x,t,oo) = e> 3 < x ' t > 3 , 

which we place on the right: 

0(A) = V(A). 5 . (4.1) 



The condition ( 2.1C ) is then automatically satisfied. We can now calculate u(x+,X-) from the pro- 
jections P(x+,X-) by using the regularity condition at A = (2.11) 



e 



20us 3 Sie -2f)us 3 =Tj ) (pyl Sl ^(p). (4.2) 



It is simple to check that the unique Hermitian projection that annihilates the vector ( 1 I is 



pi- = T+W T+1mF 
\ l+H 2 i+ImF 



The one-soliton solution 

We take one simple pole on the imaginary axis 



^(x+,x_,X) = (P ± (x + ,x^ + ^^P(x + ,x^)) (4.3) 

A — IK 

where k is real and positive. The symmetry condition ( [2.12] ) implies that 

^ X-IK ' WV X + IK 

(±±¥)(UPUi + ^— ^ UP^rf) = /(A)(P X + ^— -P). 

y X- IK ,K X + tK ' J V /K X + IK ' 

We see that 

/(A ) = (^), 

X + tK 

and 

UPU^ = P^ = 1- P. (4.4) 



The group factorisation (3T) implies that the image of P projects onto the space 

V(x+,X^) = e 2im(Kx+- rc - I x_) slFo ; 

where Vq is an arbitrary initial space. Vo must be one dimensional for there to be non-trivial solutions. 
Alternatively P 1 - is the unique projection which annihilates this space. 
We take Vq to be spanned by the vector 



\)(iQY 



I ) MY' 2 .. (4.r,i 



where Q is real. Then V{x+, xJ) is spanned by the vector 
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P^(x+,x_) = - fl _ tWs ^ W ) (4-7) 



where 

and thus the hermitian projection P (x+, x_) which annihilates V(a;+, x_) is 

n-(^) i 

This projection is also consistent with the symmetry condition ( [2.12 ) which has now been translated 
into the form (4.4). Any other choice of initial space V would not have satisfied this condition, and 
indeed it is this choice which enforces the reality of the solution. 

We have , ,„, 

/ n _rl+iWy 

#r + ,^,0) = (l-2P) = l y» 

Now 

j(j)e"^ = (1-2P) S1 (1-2P)(J 
1 -iW\ 2 i /0 



so the one-soliton solution is 



1 + iWJ 4 
1-iW 



(4.8) 



(4.9) 



1 + iW 

We note that log \ Q\ is essentially the position xq of this soliton (after suitable normalisation) at some 
fixed time, since xq is given through Q by 

sign(Q)W = e ™((«-«" 1 )t+(«+«" 1 )(^-^) ; ( 4 10 ) 

so 

log(|Q|) = -m(« + k~ 1 )xo . (4.11) 

Now Qel, and the choice Q < 0, Q > corresponds to taking solitons or anti-solitons. Note that 
we can also change the sign of k to move from a soliton to an anti-soliton. 

The position of the pole at in, is related to the rapidity 9 of the soliton by \n\ = e~ e . 
The two-soliton solution 

We take 

. i A + IKt. . . i A + IKo „ . 

^X+,X-,X) = {Pt + 7 —Pl){P^ + T — P 2 , 

A — ZKl A — ?,K2 

for «i and K2 real and positive. We can choose to write the unitary meromorphic loop corresponding 
to the pole at IK2 to the left in the product: 

Vi + T — A)(P 2 X + t r^P 2 = + T —Pz){Pt + T r^Pt ■ (4.12) 

A — IKl A — IK2 A — A — IKl 

The projections P3 and P4 are uniquely determined. We do this because we cannot immediately solve 
for P2(x + ,x_) since it is not ordered to the left, but on the other hand we can of course solve for 
P 3 (x + ,x-). 

Imposing the symmetry condition (2.12) implies 



A — ik\ A — IK2 



A + IK\ A + IK2 
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and hence, 

= /w(* + ^*)(* + i^o- 

A + ZKl A + ««2 

It follows that 

[/Pit/ 1 = Pf 1 UP 2 U^ = Pf, 
repeating the argument for the simple pole factors reversed implies 

C/P 3 C/ f = P^ UP 4 tf = Pt- 

Only the conditions for Pi and P3 are strictly needed for the calculation of P2 . Indeed we can take 
the one-soliton solution ( h? ) for Pi and P3 separately. 



1+iWi ■ 



and 



(j^-X/2 _ (^X/2 J 



where 

Wi = Q ie "»("i»+-"r 1 *-) and W 2 = Q 2 e m( ' K2X+ - K 2 lx - ) . (4.14) 



We rewrite equation (4.12) as 



A — IK2 A — IK2 



evaluating the residue at A = in 2 



U 1 (iK 2 )P 2 = P 3 U 4 (iK 2 ) 
P 2 = Ux{iK2)- l PMil^)- 

P 2 projects onto the image of Ui(iK 2 )~ 1 P-iUi{iK 2 ) and this is the image of t/i(i«2) _1 P3 since U4(ik 2 ) 
is invertible. This space is the space annihilated by P 2 . P 2 can now be computed explicitly. So we 
have 

V 2 {x +1 x-) = Ui(in 2 )~ 1 V 3 (x + ,X-), 



and it follows that 



and 



V 2 = 



l+W 1 W 2 -iY- 1 W t +iY - 1 W 2 

Y~ 1 (l-iW 1 ) 
l+W 1 W 2 +iY- 1 W 1 -iY- 1 W 2 



/ 1 1 {l+iW 1 ){l+W 1 W 2 -iY- 1 W 1 +iY- 1 W 2 ) \ 

p _ I 2 2 {l-iW^l + WiW-i+iY-^W^iY-^W^ \ (4 1*1 

2 I 1 (l-iW 1 ){l + W 1 W 2 +iY- 1 W 1 -iY- 1 W 2 ) 1 I l*- 10 J 

\ 2 (i+iH'i)(i+H'iiv 2 -jy- 1 Wi+jy- 1 iv 2 ) 2 / 

where 

K 2 + Kx 

We then find that 

C l+iY- 1 W 1 -iY~ 1 W 2 +W 1 W 2 n \ 

l-tY- 1 Wi+tY- 1 W 2 +W 1 W 2 
n i- i y-i Wl + i y- 1 w 2 +ty 1 w 2 J • 

U l+iY- 1 W 1 -iY- 1 W 2 + W 1 W 2 / 
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Now 



so we get 



4 V 1 



= ((1 - 2Pi)(l - 2P 2 ))- i s 1 (l - 2Pi)(l - 2P 2 ) 



-i/3iA 



(1 + iF" 1 ^ - iY- x W 2 + WxW 2 ) 

~ y (i - iy-Wi + iy- 1 ^ + WiW 2 ) 

If we adjust the initial subspaces Vi(0, 0) and V%(0, 0) by 
and let 



K 2 + Kl 

then we get the familiar two-soliton solution: 

_ l/3 „ _ (1 - iWi - iW 2 - XW t W 2 ) 



(1 + iWi + iW 2 - XWxW 2 ) ' 



(4.16) 

(4.17) 
(4.18) 

(4.19) 



The solution ( 4.16 ) is written in 'left-most' ordered co-ordinates, that is coordinates given by taking 
the Q's defined by the project ion f or each pole in the left-most ordered position in the factorisation. 
We note that the adjustment ( 4.17 ) to these co-ordi nates is a subtle re-normalisation of the standard 
Q's, which enter the standard form of the solution ( 4.19 ). This adjustment will be crucial for us in 
obtaining the final version of the symplectic form written in terms of the positions and momenta of 
the solitons. 

The breather solution 

We choose the poles to be at a and —a, where a is not on the imaginary axis. 



tl>(x+,x-,\) = {Pt + ^A)(P 2 ± + ^P2) 



(Pi 



X — a 
X + a 

X + a 



X + a 
X — a 



Px 



The symmetry condition ( 2.12| ) implies that 

C/PiE/t = P 3 , 

This time we take for Q G C, 



UP 2 U^ = Pa 



+ Q 



1/2 



-1 
1 



(4.20) 

(4.21) 
(4.22) 



and then, writing W(a) = e m ( ax ++ a *-) ) 

Vi(x + ,x-) = ( j) (QW»r V2 + ("/) (QW(a)) 1 / 2 , 
splitting QW(a) into a pure phase and modulus, QW(a) — UW, we have 



Vi(x + ,x-) = 



(UW) 



-1/2 



(UW) 



1/2 



((uw)- 1 ' 2 - (UW) 1 / 2 
uuw)- 1 / 2 + (UW) 1 / 2 ' ' 



(4.23) 



where, setting a — v + ik, 

U = e im( - vx + + J^+k^ = e im ^ v+ v i + k i )(*~ t o) + (-f+ ii ^" fc a ){' x ~ x o)) 
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and 



W 



e m (- fcc ++ „2 +fc a / 



e ™((- fc +^T^)( t -*o)+( fc +^^)(^- a; o)) 



Here Xo and to are fixed through the complex valued Q. We can write any Q G C — {0} as Q = 
e m(aq 1 +a q 2 ) ^ smce a nas a non-zero imaginary part, so that a and a -1 span the complex plane, and 
where q±, q2 € K. Then, on remembering that x± — t ± x, xq and to are related to q\ and (72 by 

(?i = -(%o +to), q2 = -to + x - 

Then 

Pi = 



(U-W)(1-UW) 

U(l + W 2 ) 
(-U+W)(UW+1) 
U{l+W 2 ) 



j-l+UW)(W+U) 

U(l+W 2 ) 
(U+W)(1+UW) 
U{l + W 2 ) 



and to be consistent with Eqn. (4.21), we must start with the vector 1/3(0,0): 

-1" 



^3(0,0) = Q 1/2 



+ Q 



-1/2 



1 



and we get the evolution (recall that the pole is at A = —a) 



v(x x (-{{UW)-V*-{UW)V*) 
Using the same arguments as for the two soliton case we find that, putting s = k/v, 



(4.24) 



V 2 = 



(isW+isU 2 W-U-UW 2 ) 
U(l + W 2 ) 
(UW+l)(isW+isU 2 W+U+UW 2 ) 
(~l+UW)U(l + W 2 ) 



and 



Pi = 



1 (-1 + UW)(U-W) 

2 U(l + W 2 ) 

1 {U-W)(isW+isU 2 W+U+UW 2 )(l + UW) 

2 (isW+isU 2 W-U-UW 2 )U(l + W 2 ) 



1 (W+U)(isW+isU 2 W-U-UW 2 )(-l + UW) 

2 (isW+isU 2 W+U+UW 2 )U(1+W 2 ) 

1 (1 + UW)(U+W) 

2 U(l+W 2 ) 



SO 



V(0) 



2is 



-Pi)(l 



2is 



P9 



(s 2 + l)(isW+isU 2 W+U+UW 2 ) 
' (isW +isU 2 W -U -UW 2 ){l+is) 2 



(s 2 + l)(isW+isU 2 W-U-UW 2 ) 
' (isW+isU 2 W+U+UW 2 )(l+is) 2 



The result follows: 



( 



isU- 1 +isU + W- 1 + W 



isU- 1 +isU -W- 



W 



)■ 



(4.25) 



(4.26) 



5 The symplectic form 

We have seen that the loop <fio specifies a solution to the sine Gordon equation for all space-time. We 
can consider the set of such <f>o to be the phase space of the sine Gordon system. Then a change vo in 
<j>o would represent a change in u at all values in space-time. Given two such changes vo and wq we 
should be able to calculate the canonical symplectic form (5.3), and this is what we shall do in this 
section. Here we denote the time derivative of u by u. 
Lemma 5.1 

With b defined by the factorisation (|3.l|) 



(b x b L ) v = P(5 v u)s 3 + 2m(3(S v u)(X L e 



- 1 ^ US3 [s 3 ,si]e- /3uS3 



Ae 



-f3us 3 



s 3 ,si\e 



(5.1) 
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Proof Using the results for b just prior to Eqn. ( |3.3| ) we can write 

bjr 1 = b+b- 1 - b-b- 1 = B - A 

= p(d-u + d+u)s 3 + 2m(\- 1 e )3usa s 1 e-P US3 - Xe- 0USS s x e 0US3 ) 

= f3us 3 + 2m(X- 1 e PuS3 s 1 e'^ US3 - Xe~ l3uS3 s 1 e> 3uS3 ). (5.2) 

Now vary u by a parameter we refer to as v. We obtain 

= (3(S v u)s 3 + 2m(3(S v u)(X- 1 e' 3uS3 [ S3 , Sl }e- 0US3 + Ae"^ 3 [s 3 , Sl ]e^ S3 ). □ 

The canonical symplectic form , derived from the Lagrangian density formulation of the sine-Gordon 
model, is 

02 ,00 

id = — / 8uf\8udx. 
4 J-00 

This can be written as 

01 ,00 



oj(v, w) = —— I (s v u5 w ii — S w ii5 v ii) dx (5-3) 
^ J —00 

in terms of variations v = 8 v <fi and w = 8 w <p. 
Proposition 5.2 

The form ( |5.3| ) is given by the formula 

1 f dX f°° / \ 

oj(v,w) = / —Trace / (Mr 1 ) (b^wdx (5.4) 

2m 7 7 A J-oo V ' ■» 

where the contour 7 is a small clockwise circle and a large anti-clockwise circle around the origin. Here 
'small' means that all the poles of the meromorphic function <f> lie outside the contour, and 'large' 
means that all the poles lie inside it. 
Proof 

Using Lemma 5.1, 



1 f dX f°° / 

u>(v,w) = — / —Trace / [f3(S v u)s 3 
2m J X J_ oc V 



+ 2m(3(5 v u)(X- 1 e f3uS3 [s 3 , sije"^ 3 + Xe~ 0US3 [s 3 , sije^ 3 )) tfr x w dx. (5.5) 
We calculate this integral in two parts, for the first part 



1 f d\ 

—ft l —Trace / {8 v u)s 3 (f)~ 1 w dx. 



2m .L X 

The J ^ ^ s equivalent to evaluating at A — » 00 minus evaluating at A = 0, that is 

(3 I (5„?i)Trace I s 3 <fi~ 1 w(oo) — s 3 (j)~ 1 w(0)) dx, 
J —00 ^ ' 

recall that <j>(oc) = e^ US3 and 0(0) = Ne-P US3 , where TV is some constant normalisation matrix. Then, 
we get 



/ ((5„u)Trace (s 3 )(S w u + S w u) dx 



2 r°° 

(6 v u)(5 w u) dx. (5-6) 
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The second term is 
2m(i f d\ 



Hi J T TraC<? 1" {5vU) { X ~ lel3US3 [S3 ' s ^ fjUS3 + Xe ~ f3US3 ^ ^^y^w dx, 
and integrating around the small circle in 7 gives the A coefficient of the expansion about A = of 

/oo 
( 6 V u) A - 1 Tr ace s 3 [e 0US3 s 1 e~ 0US3 , 0" 1 w] dx 
-CO 

/OO 
(^u)Trace ^[e^sie"^ 3 , (<?T ^'(O)] dx. 
-00 

To calculate this, consider the derivative of a<f>o — <f>b in the direction x_; 

A0"V- = \f3(d-u)s 3 + 2me f3uS3 s 1 e- fiuS3 - 2mf 1 s 1 (/), 
and now apply ^ at A = to get 

<TV-(0) = P{d-u)8 3 - 2m(0- 1 s 1 0'(O) - fV'f^i^O)) , 

(we use prime to denote ^) and apply -1 0_(O) = —/3(d_u)s 3 , 

-/3(<9_w)s 3 = -mO-V^O) - 0" 1 4>'<t>- 1 s 1 <j)(0)). 

Conjugating this with 0(0), and noting that 0(O)s 3 </> -1 (O) = Ns^N^ 1 , we have 

- ^(d^TVsaiV- 1 = m(</.'0- 1 (O)si -si^^-^O)) 

= m^-^O),*!]. (5.7) 

Now vary this in the direction w 

-/3(d-U w )Ns 3 N- 1 =m[w'(j>- 1 (0)-(j>'(j>- 1 w(j>- 1 (0),s 1 ] , 
and conjugate by 0(0) again get 

-0(d-u w )s 3 = m^-^'W-^-V^^W^-'^si^O)] 

= m[(0- 1 w)'(O),e' 3tlS3 s 1 e-' 3uS3 ]. (5.8) 
This gives the integral of the second term around the small circle in 7 as 

32 roc 



a 2 r°° 

— / (6 v u)(5 w (d-u))dx. 



Now wc turn to the integral around the large circle in 7. This is the A part of the Laurent expansion 
of 

/OO 
{d v u)XTr&cc s z [e- f3uS3 sie f3ua3 , cjy^w] dx 
-00 

about A = 00. We shall calculate this from the x + equation of the linear system: 

0~V+ = -f3(d+u)s 3 + 2mXe- f3uS3 s 1 e l3uS:i - 2m\<p- 1 s 1 (j>, 
taking the limit as A — > 00, and using 4>~ 1 (j) + {oo) = /3(<9+m)s 3 , 

f3(d+u)s 3 = m lim X( e - puS3 Sl e 0uaa - _1 Si0), 

A^OO 
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if we set <j>'{oo) — liniA^oo A(^(A) — 0(oo)), then 

(3(d + u)s 3 = m{4>'~ 1 4>' ' <\T X si0(oo) — _1 si0'(oo)), 
since <f>(oo) commutes with S3, we can conjugate to get 

/3(d+u)s 3 = m[0'0 _1 (oo), si], 
and differentiating along the direction w to get 

(3(d+u w )s3 — m[w'(f>~ 1 (oo) — 4>'4>~~ w<j)~~ (00), s\], 

conjugating again 

/3(d + u w )s3 — m^^w^oo) ~ (j)~ 1 cf) l (f)~ 1 w{(X)),(f)~ 1 si(j)(oo)\ 

= m[((/.- 1 w)'(cx)),e- /5tlS3 s 1 e /3uS3 ]. (5.9) 

Then the contribution to the second term for the large circle in 7 is 

52 r°° 



4 



(6 v u)(6 w (d+u)) dx 



so the total contribution from the second term is 



0i ,00 0i ,00 

(S v u)(5 w (d+u + d-u)) = — / (S v u)(S w u) dx , (5.10) 



^ ./ —00 

and the total is the required formula 

01 r°° 



cu(v,w) = — I (s v u5 w u — 6 w uS v uj dx . □ 



Now we are left with the problem of how to perform the x integral in 

1 f d\ r 



w(y,w) = /—Trace / ( b^ 1 ) dx , (5.11) 

2ttz J 1 A J_ 00 V J v 



and for this we shall need the following Lemma: 
Lemma 5.3 

-(<j>- 1 w) x + [b m b- 1 ,<j>- 1 w] = 

Proof 

Begin by varying the equation a4>o = 4>b by the parameter wq, a change in cf>o, and find 

awo = wb + 4>b w , (5.12) 
where w is the corresponding change in <j). On differentiating with respect to x we get 

a x w a = w x b + w6 K + 4> x b w + 4>b wx . 
Now substitute for wq from ( |5.12| ) and for (f> x from ( |3.2| ) 

a x a~ 1 wb + a x a~ 1 (f)b w — w x b + wb x + a x a~ 1 (f>b w — cf>b x b~ 1 b w + <frb wx , 

rearrange, 



a^a 1 w - w x - wb x b 1 = 4>{b wx b 1 - b x b 1 b w b 1 ), 
and substitute for a^a -1 from ( |3.2| ) 

(px^w + (j)b x b~ 1 (j)~ 1 w — w x — wbxb -1 = (f>(b x b -1 ) w ■ 
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Multiplying by tf> 1 on the left and rearranging gives the required result 



Proposition 5.4 



Proof 

We can write 



}) x + [b x b l ,<f>- 1 w] = {b x b- 1 ) l 



1 f d\ 

u)(v,w) — - — / —Trace b v b~ X (f)~ x w 
ziri /-, A L 



□ 



X — — oo 



1 



uj(v.w) = / — Trace 

V ' 27ri J 1 X 



(5.13) 



(5.14) 



Now integrate by parts and reorder the commutation relation using the trace property: 



u>(v, w) 



1 



dX, 



2ni L, X 



-Trace 



+ f dxbvb-^-^-W)^^ 1 ^- 1 ™])} 



Using Lemma 5.3, we can write 



lo(v, w) 



1 



dX, 



2ni L X 



-Trace 



b v b 1 <p 1 w 



+ 



dx b 



vb-^bvb- 1 ).^ , 



and we observe that the integral around 7 of b v b 1 {b x b 1 ) w vanishes since the function is analytic 
between the circles comprising 7. □ 



6 The abstract symplectic form. 

It might be thought that the form of the symplectic form we derived in the last section would be 
highly dependent on the structure of the sine-Gordon equation. However the formula ( 5.1 3| ) actually 
gives a closed 2-form associated to a group doublecross product under very general conditions, as we 
now show. 

Let X — QA4 = M.Q be a group doublecross product with an adjoint invariant inner product (, ) 
on its Lie algebra. From the factorisation a(f>o = 4>b and a change (4>q]Vq) in <fio we define derivatives 
v — (j) v — D^ . Vo ^(j> and b v = D^ OVVo ^b (in these derivatives a is kept constant). 

Our aim is to define a closed 2-form over the phase space Q with coordinate (f>o ■ To do this we first 
define another 2-form r Q on Q, for a given a £ M, as 

T a (<f>o;vo,'Wo) = ( Kb' 1 , (j>~ x w ) - ( b w b~ x , cj>~~ x v ) . 

Proposition 6.1 

dT a ((j>o;yo,VQ,w ) = ( [w q c/)q x , y <j)^ x ] , v <f>Q 1 ) - ( [0 _1 iu, <j)~ 1 y] , <\T x v ) 

- ( [b w b-\b v b- 1 ] , Kb- 1 ) . 

Proof 

First calculate the derivative 

T' a ((f)o; v , w ; yo) = ( b vy b~ x - b v b~ 1 b y b~ 1 , 4r x w ) + ( Kb' 1 , 4>~ 1( l>wv - 4'- 1 y<p~ x w ) 

— ( b wy b~ x - b w b~ 1 b y b~ 1 , <j)~ x v ) - ( b w b~ x , 4>- x 4> vy - 4>~ x y4>~ x v ) . 

By adding the other two terms and using symmetry of the double derivative, we find 

dr a {(j)o;yo,VQ,w a ) = T' a ((f) ;vo,w ;yo) - r^(0 o ; yo, w ; u ) + r^(</> ; y , v Q ; w ) 

= ( [byb^Xb- 1 ] , + ( [b^-^byb- 1 ] , ^v) 

+ ( [Kb- 1 ^- 1 ] , cj>- X y) + ( b v b~ x , [<p- x w, 4r l y] ) 
+ (b y b-\ [cp- x v,(t>- x w\) + (b w b~ x , [4>- x y,4>- x v] ) . 
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Now we try to simplify this expression, starting with the factorisation a(f>o = 4>b, and differentiating 
in the direction (<j) ; v ) to get av a = vb + <f>b v , or in a more useful form 

4>~ 1 av b~ 1 = byb^ 1 + $T x v . 

Now using this result successively, 

( [b W b-\ byb- 1 ] , ) = ( [b W b-\ byb- 1 ] , ^aVob- 1 ) - ( [b W b-\ byb- 1 ] , byb^ ) 

= ( [b w b-\ r ^yob" 1 } , </>- W _1 ) - < [b w b~\ ct>- l y] , <T W" 1 ) 

- ( [b^.byb- 1 ] , b^ 1 ) 

= ( [(j) -1 awob -1 , 4>~ 1 ay b~ 1 ] , (f)^ 1 av a b^ 1 ) 

- ( [(f>~ 1 w 1 (f>~ 1 ay b~ 1 ] , cj)^ 1 av b^ x ) 

- ( [bvb- 1 ,^] , ^avob- 1 ) - ([b^- 1 , byb- 1 }, b^ 1 ) 
= ( [wq^ 1 ,yn<t> x ] , v a <p Q 1 ) - ( [<j)- 1 w,b y b~ 1 ] , (p^avob' 1 ) 

- ([ct>- 1 w,ct>- 1 y],(t>- 1 av Q b- 1 ) - {[b w b-\r 1 y],r 1 av Q b- 1 ) 

- ( [b^.byb- 1 ] , b v b- x ) 

- ( [b^- 1 ,^] , r'v) - ( [b^- 1 ^- 1 } , b v b- x ) . 

Substitution of this into the formula above for dr a , and using adjoint invariance of the inner product, 
gives the result. □ 

Now we define a 2-form w on Q by the formula 

T 1 R 

u)((p ;vo,w ) = lim (byb- 1 , cfy^w) - ( fe^fe" 1 , cjy^v ) 

R^oc L J x=—R 

where we remember that a is a function on space-time. This is the difference in the values of r a 
between two points in space-time lying either side of the 'interesting' region, that is where the fields 
are substantially different from the vacuum. If the fields are not compactly supported, but merely 
decreasing, we take a limit 'x — > ±oo', as the points tend to regions where the fields take more vacuum- 
like values. But what does 'vacuum' mean in our group picture? A brief look at the asymptotic form 
of the soliton solutions will give a possible answer. For the sine-Gordon or principal chiral model the 
function <p(t, x) : — > GL n tends to a limit matrix valued function commuting with J and K as 
x — ► ±oo, c.f. Lemma 10.1, below. The subgroup consisting of such functions is abelian. In general, 
we assume that <j){s) tends to a limiting value in an abelian subgroup, which we call (?oo C Q- 
Proposition 6.2 

Suppose that <f)(s) in the factorisation a(s)<po = <p(s)b(s) tends to a limiting value in the abelian 
subg roup Qoc C Q as x — > oo and as x — ► — oo. Also suppose that the inner product ^, ) vanishes on 
the Lie algebra of the group M. Then dui — 0. 
Proof 

This is more or less direct from the previous proposition. The ( [wq^q 1 , yofio 1 ] > ^o^o 1 ) term in 
dr a cancels on taking differences. The term ( \byjb~ 1 , byb- 1 ] , byb- 1 ) is zero since the inner product 
vanishes on m. The term ( [cf>- 1 w 7 (f)- 1 y\ , </> -1 u) tends to zero as x — > ±oo since the commutator 
tends to zero by the abelian subgroup condition. □ 



7 The abstract higher momenta. 

Having decided that duo = under certain realistic conditions, we can try to calculate values of uj on 
certain vectors (cj>o; wo). But which vectors to use? There is a certain easy choice of vector, which will 
lead to the higher conserved momenta for sine-Gordon. 
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Proposition 7.1 

Consider c in the Lie algebra of M. which commutes with J and K 7 that is it commutes with s±. 
Then c brings about a change in (j>o through the factorisation e cr <j>o = cj)o(r)m, for some to € M.. The 
corresponding infinitesimal change wq of cf>o is given by the formula c4>o = wq + (j)od, where d is also 
in the Lie algebra of A4. Then the vector field (<fo; wo) is a Hamiltonian flow and is generated by the 
Hamiltonian 



Wo) 



2 lim 

R — >oc 



x=-R 



where d is a derivation operator acting on the Lie algebra of M., defined by Eqn. (7. 
element in the Lie algebra of M. such that dh = c. This means that 



Proof 

If we insert the formula for wq into awo 



D 



(0o;«o) 



Wo) 



(7.1) 

and h is an 

(7.2) 



we find 



b 1 aca c 



wb + i 

- bdb- 1 = ^vo + b^- 1 . 

Using our assumption that the inner product is zero on to, we obtain 

(Kb- 1 , <j)- l w) = {Kb' 1 , ^aarV) 

= ( (fibyb- 1 ^ 1 , aca -1 ) 

= ( avoipQ 1 ^ 1 — v<f)- , aca^ 1 ) 

= (uo^o 1 ) ) ~ (vcj)" 1 , aca^ 1 ) 

= ( (f)- 1 aca^ 1 <j) ~ bdb^ 1 — (j)- 1 w , <p- 1 v ) 

= ( aca -1 , ) - (bdb- 1 , ^v) - (^w , ^v) 

— (aca,- 1 , vcf)- 1 ) — (d , (j) ~ 1 vo") — ( (j)- 1 w , ) . 



(7.3) 



(7.4) 



(7.5) 



We now subtract Eqns. (7.4) and (7.5) and calculate aj(4>o;vo,wo), for some arbitrary change vq. On 
taking differences between x = ±oo, we see that the ( d, 4>q Vo ) and ( vo^q 1 , c) terms vanish as c 
and d are independent of a, so we are left with 



u((f>o]Vo,wo) = lim (<fi 1 w , 4> 1 v) — 2(v4> 1 , aca 1 ) 



x=-R 



(7.6) 



In what follows, we restrict ourselves to the case where c commutes with every a, and with the limit 



subgroup Qoo C Q. Then 



that 



zcj> - bdb- 1 ) 



as x 



6,„6 , where 



c as x — > ±oo. This means 



±oo (7r s is the projection to the Lie algebra of Q). 



Then the first term of (7.6) vanishes and we can rewrite Eqn. (7.6) as 

r i R 

u(4> ]v Ql wo) = -2 lim (v<j)~ , c) 

R— >oo L J x= — R 



(7.7) 



This might be a sufficiently simple formula to calculate oj{4>q]Vq,wq) for various vectors (4>o;vo)- 
However it is our purpose to do Hamiltonian mechanics, so we would like to answer the question 'is 
there a Hamiltonian function giving rise to the vector field (</>o;u>o) generated by c € to?'. To do this 
we would have to show that lu(4>o] vojWq) was the derivative in the direction (0q; vq) of a function of 
4>o and c. 

Suppose that there is a 1-parameter automorphism 6:KxX^X which preserves the subgroups 
Q and M, and the inner product on the Lie algebra of X. Then there is a derivation^ $ on the Lie 
algebra defined by 

%)=e"(0,e;0,y;l,0) (7.8) 
(e is the group identity), and this preserves the inner product, that is (fly, z) + (y, $z) = 0. 

1 Loosely speaking, in the applications, i? can be thought of as ■& = X- 



d 

' dX 
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If there is an ft in the Lie algebra of M. so that dh = c, then 



uj{(/)q;v ,w ) 



-2 lim 

R^oo 



(vc/)- 1 , tih) 



x=-R 



2 lim 

R^oc 



x=-R 



(7.9) 



For our purposes we wish to to swap the order of $ and Dr<j, - Vo \ i n this expression. To do this we 
proceed cautiously and apply i? to vcf)^ 1 to get ^(v)^ 1 — t>0 -1 ^^ -1 , where $(«) = 6"(0, (j>; 0, u; 1, 0) 
and <fi§ = G'(0, 0; 1,0). From the last formula and the symmetry of double derivatives we see that 



tf^" 1 ) 
(tf^" 1 ), ft) 



D 



(</>o;«o) 



[v0 1 , 0,90 X ] , 



D(4, q - Vq ){M 1 ,h) - ([v(j> l ,<j><s<f> x ] , ft) . 



If $ preserves the abelian subgroup Goo then the Lie bracket \4>^4> 1 , vcf> 1 tends to zero as a; — > ±oo, 
and we can write Eqn. (7.9) 



w(0o;v o ,u> o ) 



( 0tf0 1 , ft) 



or in other words the vector field (0o! w o) generated by an element c G m is a Hamiltonian flow, with 
Hamiltonian 

i-i 



2 lim 

i?->oo 



R 



□ 



x=-R 



Now that we have expressions for Hamiltonian functions f c we should show that these functions 
Poisson commute. 
Proposition 7.2 

The Hamiltonians {/ c } Poisson commute. 
Proof 

Suppose that the vector (cj)Q;vo) is generated by another element c in the Lie algebra of M which 
commutes with si. Then c0o = vq + <f>od, where d is also in the Lie algebra of Ai. Eqn. (7.7) applied 
for an arbitrary we can insert this particular vq to find, 



w{<j> ;v ,wo) 



2 lim 

R^oo 



(v^ 1 , c) 



x=-R 



-2 lim 

R^oo 







z=-.R 



This means that the function f c Poisson commutes with /g, and therefore the {f c } are in involution. 

□ 

As we shall see, in the case of the sine-Gordon equation, we can take c = J or c = K to obtain 
the Hamiltonians representing the total energy and momentum in light-cone co-ordinates. The higher 
momenta are obtained by taking the infinite number of other possible c's. Since the energy is the 
Hamiltonian generator of time translations, and the momentum is the generator of space translations, 
all the higher momenta are conserved by these translations. The 1-parameter flow corresponding to 
the Hamiltonian f c on the phase space is given by the factorisation (r, <fio) t— > (f>o(r), where 4>o{r) G Q 
is the solution to the factorisation problem e rc 0o = 0o( r )^ for some d £ M. 



8 The higher momenta for sine- Gordon. 

Here we shall specialise the results of the last section to the solitons in the sine-Gordon model. The 
1-parameter automorphism of the loop group X is given by 0(s,p)(A) = p(s\) for A G C*, which is 
actually the Lorentz boost for the system, giving $ — A^. The inner product is 

1 f d\ 

(y,z) = — Trace / —y{X)z{\), 
zm J A 



15 



which is 9 invariant. We can choose h(X) — siA™, giving c(A) = nsi\ n . The meromorphic loops for 
sine-Gordon split into two cases, solitons and breathers, and we shall calculate the higher momenta 
for both cases. 

The loop <j> for a single soliton is <\> — Nipe^ US3 , where N commutes with si, and 

A + in 



for k real. Then 



i>= P 



-2inPsi \ r - 



A 



- — Trace / 
2m J 1 (X — lK){X + iK) 



IK 



dX 



-(^"(l + ^l)"- 1 ) Tracc(Psi) , 



which is zero if n is even, and if n is odd we calculate 



fc = -4/C n i(-l)(' 



-l)/2 



Trace(Psi) 



= -4«; n i(-l)( T ' 



-l)/2 



4 14 



and inserting the x dependence of fx shows that 

The calculation for the breather is more complicated. In this case 
by ( 4.20| ). We write ip = tpifa, where 

A - 



P 



Pi 



X 



and tf>2 



P, A 



Pi 



= tpe l3uS3 , where tp is given 
A + a> 



A + a 



Now tfO)^ 1 = tfOOV'i" 1 + ^(^2)^2 Vr 1 , and we can calculate 



(0(V>i)C\fc) = ^ Trace/ 



-(a — a)P\Si X r 
(X-a)(X-a) 



dX = -(a n ~ a") Trace(PiSi) . 



The term (ipi^ifc)^ bas the same limit as x — > ±00 as the simpler term ('&(tp2)tp2 

because ipi tends to a limiting loop which commutes with s\. Now we can write 



fc 



{${<t>)<t>-\h) 



X— — OO 



-2{a n - a n ) Trace(PiSi) + (-l)"Trace(P 2 si; 



x — — 00 



For the breather case we can assume that a lies in the upper half plane, in which case the equation 
above gives zero for n even, and for n odd, 

f e - 2i(a n - a n ) . 

The integer n of the n th conserved charge is the Lorentz spin. The fact that these charges were 
zero for n even was observed beforehand in in the more general context of the afHne Toda field 
theories, since the elements of the principal Heisenberg subalgebra of the Lie algebra of M, only 
has odd principal grades. In || it was shown that elements of the principal Heisenberg subalgebra 
generated the higher conserved charges, with the derivation operator in the principle grade acting as 
the Lorentz boost, indeed the principal Heisenberg subalgebra is nothing more than the subalgebra 
of the Lie algebra of M. which commutes with J and K. It is then an immediate consequence that 
the Lorentz spins of the charges, which are measured by the values of the principal grade, are only 
non-zero for the exponents of the affine algebra to. For m = su(2), these are the odd integers. The 
argument presented here is a concrete verification of this. 

We would also like to bring to the reader's attention the expressions given by the central parts of 
<\T X J(f>, and 

< J^O)^ 1 > ,and < K^i^cj)- 1 >, 

respectively, discussed in our previous paper [0, which were identified as the energy and momentum 
densities in light-cone co-ordinates integrated up to a point x. It is evident that these are the same 
as the Hamiltonian expressions (7.1), for c = J and c = K, provided we evaluate the differences as 
x —> ±00 of the integrated expressions. This is because the Hamiltonian expressions only give the 
total quantities. 
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9 The splitting of the symplectic form. 



We must now address the practical concerns of calculating the symplectic form. First we deal with 
the normalisation, that is (f> = N^pe^ US;s , where tp = 1 at A = oo, and where the constant matrix N 
(which commutes with si) may be needed to satisfy the symmetry condition. We wish to rewrite the 
expression for the symplectic form in terms of ip, which is easier to deal with. After doing this we 
see that the form splits into two parts, one due to 'self interactions and the other due to 'mutual' 
interactions. 

It is not too difficult to see that N cancels immediately from the expression for the symplectic 
form, as it cancels from ^r x (\) Wl and does not affect b at all. We shall therefore continue assuming 
that N = 1. 
Proposition 9.1 

Consider the pole at a, we write ip = £%, where £ is of the form ( 2.1 3| ) having a pole at a, and x 
is regular at a. Perform the usual factorisation (3.1) on £ only, that is a^o = C^i this defines d € M. 
Then the a contribution to the symplectic form (5.13) is 

r i 00 

u(v,w)\ a = (dQ rw ^-r'c^r 1 ^) - <x»x _1 ,r 1 c«;> Qa ■ (9.1) 

The notation <, > a ,a means take the contributions to the integral in the inner product only from the 
poles at a and a. 

We can interpret the first term here as a self interaction term of the a pole with itself, and the 
second term as the mutual interaction of the a pole with the other poles. The total form is evidently 
the sum of such terms for each of the poles a. 
Proof 

First note that 

and 

b v b~ l = Puovbssb- 1 + e- /3uS3 (^- 1 a^o^o" 1 a"V-V'" 1 ^)e /3 " S3 - 0u v s 3 , (9.2) 

Now loo k at the contribution to (b v b~ l , t/) -1 ^), which makes up the form ( |5.13 ), from the Pu^ybs^b -1 
term in ( |9.2|) , which is 

(/3uo„6s36 _1 , <jr x <j) w ) = (3u v{s3, b' 1 4>~ x 4> w b) = Puo v (ss, <t>o 1 (j>ow) - /3u 0v (s 3 , b~ x b u ) . 

Here the second term is zero as the inner product vanishes on analytic functions, and the first term 
cancels on taking the difference of its values between x = R and x = —R. The contribution to the 
inner product (fe^fe -1 , from the (3u v s 3 term of ( |9.2| ) tends to zero as x — > ±00, and we are left 

with the contribution 

which simplifies to 

and the part containing f3u w s 3 vanishes as it is the inner product of two analytic functions, giving the 
simpler result 

Here we have used the factorisation o/0o = ip c - As cipQ 1 ipo v c~ 1 — ip~ 1 'tp v is analytic, we only have to 
sum over the residues of ip~ 1 ip w to calculate the inner product. 

Let us look at the contribution from the poles at a (and a) to the inner product. We suppose that 
%jj = (x, where Q has a pole at a, and x is regular at a. Then we can write 

^~ lr ipw = X~ 1 Xw + X^ 1 C 1 CwX ■ 
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If we also write a£o = then c = \ "XO) and the a and a contribution to the inner product is 

+ (dQ 1 ^ 1 - C 1 Cv,C\ w ) a 5 (9.3) 



The first term in (9.3) is 



(XQvX Q l ,d 1 C 1 ^d) a & = (xovXo 1 ^o 1 Cow) a & - {xovXo 1 ^ 1 ^) tt , a > 

the last term of this vanishes as there is no pole at a or a, and the first term vanishes on taking the 
difference between the two asymptotic values of x. We can now write the a-contribution to the inner 
product as 

(dCo'Covd-'-C'UC 1 ^)^ - (xvx^KC 1 ^)^, 

which can be rewritten as 

as C _1 C«i only has poles at a and a. □ 

10 The mutual interaction term. 

The limiting behaviour of the meromorphic loops. 
Lemma 10.1 

The projections Pi, making up the product of meromorphic loops of the form (2.13) for </>(A), 
commute with si in the limits x —> ±oo. 
Proof 

From the equations of motion of the linear system 

4>~ 1 d + (i> = A-<l>- 1 J<i>, and _1 d_<^ = B — $~ X K<$>. 

in the limits x — ► ±oo, (for finite t), the left-hand sides of these equations are zero, and A — > J, 
B — » K, so we conclude that <^(A, ±oo) commutes with s\. We add a soliton to the system by 



multiplying a previous by a meromorphic unitary loop on the right of the form (2. IS), and adjusting 
the normalisation. The space-time behaviour of the previous <f> is unchanged. We also find that the 
limit of the resultant <f> commutes with s±, so we see that the limit of the meromorphic loop we have 
added also commutes with si. By taking residues therefore any of the Pi in the products of the form 



(2.13) making up <f> also commute with s±. □ 

For the following, recall that all matrices commuting with s± commute amongst themselves. 
Lemma 10.2 

We factor ip into the form ip — C/3X/3; where (,p has only a pole at ft, and xp is regular at (3. £p is 



of the form (2.13), where we denote the relevant projection by Pp(x). Then for any two poles /?i and 
(3 2 in ip 

Trace(P ft (oo)P fe (oo) - P 01 (-oo)P A (-oo)) = . 

Proof 



This result is almost immediate from the form (4.7) of the projection for the one-soliton solution. 



In the limit (x — > — oo, Im (f3) > 0) or (x — > oo, Im (0) < 0), 
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and for (x — ► oo, Im (8) > 0) or (x — > — oo, Im (/?) < 0), 



P/3 



1 / 1 



-1 1 



A simple calculation establishes the validity of the statement. □ 
Proposition 10.3 

The mutual interaction term, the second term of ( |9.l| ), vanishes, i.e. 



{XavXa >Ca Caw) a a 



= 



x — — oo 



Proof 

Given a unitary meromorphic loop ip with ^>(oo) = 1, we can factor it into 

i> = (0x0 , 



(10.1) 



where C/3 has only a pole at /?, and xp is regular at B, but will contain a pole at a, c.f. equation ( |9.l| ), 
unless we choose 8 = a. In the generic case where poles are simple, 



A-/3 
A-/3 



Pfl 



c,- 1 = p, 



A — j9 



P 



/3 > 



0-£ 
X-B 



PpPpv 



; A^ 

KPv -Pv)+Pvf3-Pj 

(X-3)(X-0) 



Pfl 



(10.2) 



where B v and P^ are the changes in 8 and Pp, respectively, due to the vector Vo- In deriving this 
formul a we must remember that P^Pv = P V P, which comes from differentiating the equation P 2 = P. 
From ( 10.1 ), we compute 

= Xp X Cfi l CpvXp + Xp X Xpv , (10-3) 

and noting that ip~ 1 ip v is a meromorphic function, zero at infinity, and has simple poles at 8 and /3, 
we can write if)~~ 1 ijj v in terms of its simple poles as 



= \ - / res^" 1 ^) res^(^ V«) 
V Vv x-B 



X-B 



From ( 10.5 ) this is equal to 



! _ >p f XpiP) 1 res^ 1 C/3v) xMfi) j X/3(/3) 1 resg(C /3 1 C,f3v) Xff(P) \ 

V V>v- X _Q + X _g ) 



(10.4) 



We re- arrange formula ( 10.3 ): 



(10.5) 



and c alcula te Xa vXn 1 ( a )i which we will eventually insert into < XavXa ■> (a 1 (aw >a,&- Insert the 
form (|10.4[) in to (|l0.5|) , and consider the contribution to the sum from the 8 ^ a poles. 
From Eqn. ( 10.2 ), and noting that Pp v — ► 0, as x — > ±oo, res^^J Xpv) — * /3i>P/3(±oo), and 
res^(C^" 1 C/3«) ~~ * ~~ A;P,s(±oo), as a; — > ±oo, and since x a (a)x/3(/3) _1 , and Xa( a )X/3(/^) _1 tend to matri- 
ces which commute with si, c.f. Lemma 10.1, and therefore commute respectively with res^C^ C/3„), 
and res ( g(C ( ^ 1 ('/3 1 ,) after the limits x —> ±oo, we get the contribution to XouXqH ) °f 



A, 



a — 8 a — 8 



Pp(±oc) 



(10.6) 
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The contribution to XavXa 1 ^) from the a pole is zero, and for the a pole we get the limit as A — > a 

of 

XaWXa(a) 1 res a (C a 1 C^) X«(a)Xa(A) 1 -res a (C a 1 Ca^) 

A — a 

This is 

[Xa(A) XQ (a)~ 1 , res^C^a,)] + 0(A - a) , 

which tends to zero on putting A = a and taking limits as x — > ±00, since after these limits, both 
terms in the commutator separately commute with s l5 c.f. Lemma 10.1, and so commute with each 
other. Thus XavXc^i 01 ) ^ s equal to formula ( 10. 6| ), at the limits when x ±00. In the formula 
( 10. 2| ) for Ca 1 CawW we note that Pp v — > 0, as x — > ±00, so that only the last term of ( 10. 2| ) 
contributes. Since XauXa^A) is regular at a, we see that the a pole contribution to the inner product 
(XavXa 1 ) Ca 1 C<m«) q a is made up of sums of terms proportional to 

trace(P ( a(oo)P ct (co) — P^(— oo)P Q (— 00)) 

for P ^ a. From Lemma 10.2, these terms are all zero. Likewise the contribution from the a pole 
gives zero, so we conclude that the mutual interaction term (xavXa 1 ^ Ca 1 (aw) a - contributes nothing 
to the symplectic form. □ 



11 The self interaction term. 

Proposition 11.1 

Consider the self-interaction term, < dCo~ 1 Co«^~ 1 ~ (~ 1 Cv, C _1 Cuj >> the first term of ( |9.lD , then 
we have 

(dCo-'Co^d- 1 -C^C 1 ^) = ±2Re UvQw ~ awQv , (11.1) 
(where + corresponds to Ima > and — to Ima < 0), or in wedge notation, 

= ±2Re daAd Q. ( 11. 2) 

aQ 

Here Q is the coefficient specifying the initial projection of the left-most ordered meromorphic loop 
in the product of loops. For the two-soliton case, comparing with Eqn. (4.14), this Q is either Q\ 
or Q2 depending on whether the pole at in\ or ik% is considered. The total form is the sum of the 
contributions of this type. 
Proof 
We take 

c = + , 

A — a 

which gives, as before, 

CX = ^P,P + ^PP, + M°« -*0 +*»°-^* p , ( n.3) 
A — a A — a (A — a) (A — a) 

and for convenience we shall write n = X(a v — a v ) + a v a — a v a. We also need to look at the value of 
d = C _1 <o at a and a. We recall that C" 1 = (P x + $E§P), and then 

d(X) = (P*- + ^P)a(X)(P x + ^P ) . (11.4) 

A — a A — a 



Now, from (11.3), we find that the first entry in the inner product to be calculated is 
H(X) = dC^Covd- 1 - C% 



a — a 



, (dPovPod- 1 P V P) + ^(dPoPovd- 1 - PP V ) + "[ dP0d ~\ P l 

A — a A — a (A — a)(A — a) 
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The contribution to the total inner product from the a and a poles, noting that P 2 = P, using 
cyclicity of the trace, and P ± P W = P W P, is: 

/ ,a w , N a w , ^ \ , N / ,PH(a)P x P ± H(a)P s \ 
Trace(P(— H(a) - -^H(a))P) + (a-a)Trace(( ^ + )P W ) , (11.5) 



a a / V a a 



and we calculate, by substituting (11.4) for d(X), and writing a for a(A), 



p „ mp , ^( PaPovPoa-'P , PoPpP^^^ 

PH{\)P = (a-a)( + - j + n— — — 

\ A — a A — a / {A — a){A — a) 

Dff/UDi , _./(A-a)Pafl h ,fl,a- 1 P ± , PaP P 0v a- 1 P ± - PP V \ PaP^P 1 ^ 
PH(X)P = (a-a)( ^— ^ + — J + " (A _ a)2 

Dicf , UD , _>(P ± aP 0v P a~ 1 P-P v P {X-a)P ± aP P 0v a- 1 P\ P^aP^P 
P H{X)P = (a-a)( ^ + ^-^ j + » (A _ a)2 ■ 

We know that P is the orthogonal projection to a(a)imp), which gives rise to the equations 

Pa(a)P = a(a)P Q and P a{a)~ l P = a(a) _1 P . (11.6) 
Using these equations ( [11 .6 ) , we find 



PH(a)P = PaP 0v P a- L P - (a - a)PaP Po v a- L a x a- L P 



a v (Pa x P a- 1 P - PaPaa^axa^ 1 P) 



If we take the equation Pa(a)Ppa(a) 1 P = P, derived from (11.6), and differentiate it, we find that 

P v aP Q a- l P + PaPoya^P + PaP^Py + a^PaxPoa^P - PaP^ax^ 1 P) = P v , 

and using this we can write 

Trace(PPT(a)P) = - Trace (aP a _1 P^ + (a - a)PaP P 0v a- 1 axa- 1 P). 

We also find that 

PiJ(a)P ± = aP P 0v a~ l P x - PP V + -^^aPoa^P 1 - 

a — a 

Trace(PPT(a)P- L P u ,) = Trace (PaP P vP^a' 1 P w ~ PP V P W + -^^aP aT l P W P) . 



The combination from the X — a terms of ( 11.5 ) is, writing a for a(a), 

Ti&ce{a w PH(a)P+(a-a)PH(a)P ± P w )/a = a _1 x 

Trace(a„aPoa- 1 P w P-a li ,aPoa- 1 P u P+(a-a)(PaP Po u P ± a- 1 P w +a w PaP Po„P ± (a- 1 )AP-PP u P ii -))- 



By differentiating the equation P^~a 1 (a)P = 0, derived from ( ll.q ), we obtain 



Pou.a-'P + a w P i (a- 1 ) A P + P - L o- 1 P„ = , 



so 



Trace(a w PH(a)P + (a - a)PH(a)P ± P w ) /a = 

TTa.ce(a v aP a- 1 P w P - a^Poa^PyP + (a - a)(P P Qv P 0w - PP v P w ))/a . 
The term PP V P W tends to zero as x — > ±oo, and PqPqvPqw vanishes on taking the difference between 



two values of x, leaving the X = a terms from (11.5) as 



i [Trace(a !i ,PP(a)P + {a - a)PH(a)P ± P w )] ~ _ oo = 
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^■[Tmce(a v aP a- 1 P w P - a^Poa^P.P)]^ ^ . (11.7) 
Let P be the orthogonal projection to the complex vector (1,/i). Then 

P=^—( l M and P,P = ^ (-* 



l + Hii) ' (1 + MM) 2 V 1 M 

Since a(A) = e - Jx +- Kx - j it can be written in the form 

1 ( r + 1/r r — 1/r 



v 7 2 \ r - 1/r r + l/r J ' 

where r = exp ( — imx(a — a~ 1 )/2^. Here we have set t — 0, although the limits derived below are 
equally valid for any finite t. Then it is possible to calculate 

TraceKaPoa-P.P) = — , ^j^P^fc^ w , ^ a, r ' 



(1 + |mo| 2 )((1 - Mo)(l - Mo) + M 4 (l + Mo)(l + Mo)) ° r 

qig-i)(H 4 -i) 

(1 + |moI 2 )((1 - Mo)(l - Mo) + k| 4 (l + Mo)(l + Mo)) 



As a v r w /r = — imx(l + a 2 )a v a w /2, the first term vanishes on subtraction of the term with v and 
w reversed. If Im(a) > 0, so that the lower limit corresponds to \r\ — 0, equation (|ll.7|) becomes 

-[Trace(at,aP a Pu>P - a w aP Q a P V P)\ , r , =0 = , 2 tt ■ (H-°) 



The overall sign of (11.5) is reversed if Im(a) < 0. If we use the coordinate Q, c.f. Eqn (4.5), where 

1-iQ 



Mo 



1 + iQ 



we can rewrite Eqn. (11.8) as 



a v Q 

a 



— rTrace(a«aPoa 1 P w P — a w aPoa 1 PvP)]T'i—n 



We also calculate the A = a parts of (11.5), and find the total contribution to the symplectic form 

(d( Q L ( 0v d L -C Cv,C Cw) = ±2Re — , 

where + corresponds to Im(a) > and — to Im(a) < 0. □ 



12 Final statement of result. 

The two-soliton case. By combining Propositions 9.1, 10.3 and 11.1, and by considering the 
contributions from the poles at in\, and ZK2, we have derived the symplectic form u) for the two- 
soliton solution 

lo dni A dQ[ dn^ A dQ' 2 
2 KiQ[ K2Q2 

Throughout this section, we must take the real parts of the expressions giving the forms. Also, for 
the solitons, without loss of gene rality we can take k > 0, and therefore take the plus sign of ( 11. 2j ). 



If k < 0, the minus sig n in (11.2 ) is to be expected because of the relation of the position xq of the 



soliton to Q, see Eqn. ( |4.11|) . If we wished to write the form in terms of the rapidity 9 and position 
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of the soliton, we must substitute \k\ = e 9 , and therefore pick up a minus sign if k < 0, this cancels 



with the minus appearing in (11.2), so in both cases we would have an overall plus sign. 

(12.1) 



Here the explicit two-soliton solution in these co-ordinates is written as 

_ (1 + iY-iC&Wi - iY- x Q' 2 W 2 + Q'iQ'zWxWz 



(1 - iY-iCfiWi + iY-^Q'JNi + Qi<%WiW2 
with 

_ K 2 - Ki 



compare with Eqn. (4.16), and following. Also, Wi = e m ( KiX + K t x ~\ and Ki is related to the rapidity 
9i of the soliton by e~ tii = \kA. 



If we write (12.1) in terms of the co-ordinates Qi,Q 2 which are more familiar to us, so that the 



two-soliton solution is 

_^ = (1 + iQxWi - iQ2W 2 + Y 2 Q 1 Q 2 W 1 W 2 ) 
6 ~ {l-iQ 1 W 1 +iQ 2 W 2 + Y 2 Q 1 Q 2 W 1 W 2 )' 

Qi = F- x Qi, Q 2 = F- X Q^. 



by defining 



Then the symplectic form is 

w _ d^AdiYQx) dn 2 Ad(YQ 2 ) 
2 = KiYQi + k 2 YQ 2 

dni A dQi dn 2 A dQ 2 dn\ A dY dn 2 A dY 
KlQl k 2 Q 2 KlY n 2 Y 

However, it is easy to show that 

dK X A dY dn 2 A dY 2 



(12.2) 



K\ K 2 («1 + K 2 ) 

so we can rewrite the form as 

ui dni Ad(Y 2 Q 1 ) dn 2 AdQ 



dni A dn 2 



2 kiF 2 Qi k 2 Q 2 

or 

w d«i A dQi dn 2 A d(Y 2 Q 2 ) 

2 = KiQi + At2^ 2 (92 ' 

these last two expressions are ones familiar from Babelon and Bernard [jj], where the form is rep- 
resented in terms of "in" and "out" co-ordinates, respectively. The multiplication of one of the Q's 
by Y 2 , is precisely the shift in one of the co-ordinates of the solitons by the standard time delay, 
proportional to log(F 2 ), for sine-Gordon. 

The n-soliton case. The form is evidently the sum of single soliton contributions given by the left- 
most ordered Q's, denoted QL However, it must be remembered that these Q\ are not the co-ordinates 



which appear in the standard n-soliton solution, that is the n-soliton version of ( 4.19 ), or the n-soliton 
solutions given in |H, (in they are denoted by Of). We denote these co-ordinates Qi. The two sets 
of co-ordinates are related by 

q: q\\ 



Exactly as we have seen for the two-soliton case above, we can then make contact with the diagonal 
result of 0], where the form is given in terms of 'in' or 'out' co-ordinates. 
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The breather case. Begin with the factorisation (4.20), where the left ordered projections are Pi 
and P3. Without loss of generality we may assume that the imaginary part of a is positive. Then the 
contribution to the symplectic form from the a pole is 

da A dQi 



To find the contribution from the —a pole we use equation (4.21) and (4.24) to find Q3 = l/Qi, and 



we can now give the total symplectic form for a breather, remembering to take the real part of this, 
as 

da A dQi 



= 2- 



aQi 



The Q used here is the one appearing in the left-most ordered factor, the factor with a pole at a, in a 
product of loops representing many breathers and solitons. It is crucial to remember that this Q will 



be shifted, as compared with the Q in the single breather solution (4.26), by factors which arise when 
the other breathers and solitons to the right of the left-most ordered factors are taken into account. 
Compare with the n-soliton case. 
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